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Abstract We demonstrate that magnetic 3 d impuri¬ 
ties with S = 3/2 spins and no orbital degree of free¬ 
dom induce changes of spin-orbital order in a 4 d 4 Mott 
insulator with S = 1 spins. Impurities act either as 
spin defects which decouple from the surrounding ions, 
or trigger orbital polarons along 3d-Ad bonds. The Ad- 
Ad superexchange in the host Jhost competes with 3d- 
4d superexchange Ji mp — it depends on which orbital 
is doubly occupied. The spin-orbital order within the 
host is totally modified at doping x = 1/4. Our find¬ 
ings provide new perspective for future theoretical and 
experimental studies of doped transition-metal oxides. 
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Fig. 1 Left — spin-orbital C-AF/G-AO order found in Mott 
insulators with d 2 (vanadates) or dA (ruthenates) ionic con¬ 
figurations. Spins 5 = 1 are shown by red (black) arrows and 
doubly occupied f 2 g orbitals (doublons) are shown by green 
symbols for a and c orbitals. Right — the orbital dilution oc¬ 
curs after doping by a 3 d 3 ion with 5 = 3/2 and no orbital 
degree of freedom. Here a doublon site is replaced by Mn 4 + 
or Cr 3 + impurity coupled to the host by 3d-4d bonds. 


1 Orbital dilution 

Entangled spin-orbital superexchange interactions lead 
to several surprises in transition-metal oxides [jQ. In 
Mott insulators these interactions are modified by dop¬ 
ing which may generate novel phases — they emerge 
from the interplay of complex spin-orbital-charge cou¬ 
plings. For instance, hole doping of a ferromagnetic 
(FM) system with d 1 ionic configurations removes lo¬ 
cally orbital degrees of freedom and generates stripe 
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phases with orbital polarons [2]. We show below that 
similar polarons emerge also in doped spin-orbital sys¬ 
tems. 

The well known example are hole doped Lai-^Sr^MnOa 
manganites with colossal magnetoresistance [Sj ■ At low 
doping orbital polarons emerge in an antiferromagnetic 
(AF) system by double-exchange mechanism |3], while 
at higher doping spin order changes globally to a FM 
metal coexisting with an e g orbital liquid phase [5]. But 
a different scenario is also possible — frustrated spin- 
orbital interactions may lead in some cases to the col¬ 
lapse of any long range order and to a spin-orbital liq¬ 
uid suggested for LiNiC >2 [6]. However, spin and orbital 
energy scales are here quite different and the reasons 
behind the absence of magnetic long range order are 
indeed more subtle and not yet fully understood [7 . 

A rather unique example of a spin-orbital system 
are perovskite vanadates, where interesting competi¬ 
tion between two types of spin-orbital order was ob- 
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served [ 8 ]. The theoretical spin-orbital model describes 
an interplay between 5 = 1 spins and r = 1/2 orbital 
doublet {yz, zx} active along the c axis [5], while xy 
orbitals are filled by one electron each. In addition, G- 
type AF (G-AF) order in YVO 3 is fragile and switches 
to C-type AF (G-AF) one in Y|_. r Ca x V 03 [TO], with 
staggered lines of one-dimensional (ID) FM order (Fig. 
1), accompanied by G-type alternating orbital (G-AO) 
order. Already at low doping x ~ 0.02 finite spectral 
weight is generated within the Mott-Hubbard gap due 
to charge defects away from the VC >6 octahedra HU. 

Here we consider doping by static charge defects, 
e.g. Mn 4 + or Cr 3+ ions, within the ab planes having 
columnar G-AF order of S = 1 spins at ruthenium or 
vanadium ions in the host which leads to orbital dilu¬ 
tion [12], see Fig. 1. In contrast to hole doped man- 
ganites, holes are here immobile and disturb t 2 g orbital 
order. It has been shown that: (*) dilute Cr doping for 
Ru reduces the temperature of the orthorhombic dis¬ 
tortion and induces FM behavior in Ca 2 Rui_ a ,Cr a; 04 
(with 0 < x < 0.13) [13]; (ii) Mn-substituted single 
crystals of Sr 3 Ru 2 _ x Mn ;c 07 reveal an unusual E -type 
AF structure at x = 0.16 mi which is again triggered 
by double exchange. These findings motivate the theo¬ 
retical search for the consequences of orbital dilution. 


2 Isolated impurities 




Fig. 2 Phase diagram for the 3d impurity in the ab plane 
of 4 d host with C-AF /G-AO order and impurity replacing a 
orbital doublon. Spin order is switched on the central vertical 
line in phase AFal (dashed boxes), while in FSal phase the 
impurity spin is frustrated and unable to alter C-AF pattern 
in the host; at large r?i mp FM order is found in FMal phase. 


For CaoRuCH host we use the spin-orbital model de¬ 
rived for t\ g ions with low (S = 1) spins : T5] . This 
model uses t 2g doublons as orbital degrees of freedom 
and is isomorphic to the vanadate d 2 model |SJ, with the 
doublons transforming into empty orbitals (occupied by 
two holes). We shall label t 2g orbitals by index 7 when 
a given orbital is inactive along a direction 7 £ {a, b , c}: 

I a) = | yz), 16) = | xz) , ]c) = | xy ). (1) 

We consider a two-dimensional (2D) square lattice with 
transition-metal ions connected via oxygen orbitals as 
in an ab RuC >2 plane of Ca 2 Ru 04 (S 1 RUO 3 ). In this 
case |a) (| b)) orbitals are active along the b ( a ) axis, 
while |c) orbitals are active along both axes, a and b. 
The superexchange in the host for the bonds (ij) along 
the 7 £ {a, 6 } axis, 

n id - Ad = j^ ost Y {JijHSi-Sj + i)Jij + k\] ] Y ( 2 ) 

(b')ll7 

is given by Jh os t and depends on orbital operators, j[p 
and Kp } jig. The above form is generic [I] and the 
interactions depend on the intraorbital Coulomb U 2 el¬ 
ement and Hund’s exchange J 2 in the host. 


We introduce two parameters to characterize the in¬ 
teractions along the impurity-host 3d-4d bonds i f?; : 

1 - — - ,J 1 

Jimp — ! ?7imp — ^ • 1A ) 

Here the charge excitations 3df 4d 4 => 3dfM 3 determine 
the d 3 -d 4 superexchange Ql and involve the energy 

A = I e + 3 ( 1/4 — U 2 ) — 4(Ji — J 2 ). (4) 


It depends on the onsite Coulomb interactions {U m } 
(rn = 1 stands for the impurity d 3 ion), on Hund’s 
exchange {J m }, and on the ionic energy I e . For Mn or 
Cr impurities in ruthenates A > 0. 

With the parametrization introduced above, the dom¬ 
inant term in the impurity-host Hamiltonian for the 
impurity spin S; interacting with the neighboring host 
spins {Sj} for nearest neighbors j £ J\f(i) can be writ¬ 
ten in a rather compact form as follows 

%d-4 d (i) - Y { J s(Dp ] )( Si-Sj) + E d dW} , (5) 

7,ieA f(i) 

where the spin couplings Js(Z)j 7 ^) depend on orbital 

(doublon) configuration, Dp' 1 is the doublon projection 
operator at site j and the doublon energy Ejj depends 
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on ??imp Eq. (3). It can be shown [12] that the domi¬ 
nant energy scale is EP d , so for a single 3d-Ad bond the 
doublon does not occupy the inactive ( 7 ) orbital and 
spins couple with Js(Dj y ' = 0) which can be either AF 
if t/iinp < 0.43 or FM if 7 y imp > 0.43. The change of 
sign at rfi — 0.43 is reminiscent of that found in the 
Kugel-Khomskii model and as there [TB] could lead to 
exotic spin phases. 

It has been found that a single d 3 impurity at site i 
modifies the spin-orbital order at its nearest neighbors 
j £ A/"(z), while second nearest neighbors were assumed 
[121 to follow the G-AF/G-AO order in the host (Fig. 1). 
Here we release this constraint and consider classically 
the impurity with its first, second and third nearest 
neighbors. One finds that the phase diagram is almost 
unchanged by increasing cluster size from | 12 | for the 
doping at c doublon and doublon orbitals at nearest 
neighbor sites change from inactive to active ones with 
increasing J lmp (not shown), the latter similar to orbital 
polarons in manganites Ba¬ 
in contrast, for a doublon doping the impurity flips 
a to & orbitals at second neighbors along the vertical 
line and reverses four spins along the b axis (replac¬ 
ing a ID FM order along if) in phase AFal for small 
r/i m p, see Fig. 2. This shows that the modification of 
spin-orbital order induced by 3d-Ad bonds may be long 
range. Surprisingly, the host G-AF/G-AO order recov¬ 
ers for larger 77 i mp ~ r/ mp in frustrated spin FSal phase 
when the AF 3d-Ad coupling weakens. Here the quan¬ 
tum fluctuations support the impurity spin following 
G-AF spin order. Finally, when the impurity-host su¬ 
perexchange becomes strong FM, the entire impurity 
cluster polarizes ferromagnetically in FMal phase. 


3 Phase diagram at x — 1/4 doping 

As for a single impurity, the 3d-Ad bonds strongly influ¬ 
ence spin-orbital order at finite doping except for a rela¬ 
tively narrow window of Hund’s exchange ? 7 im P — ^fmp- 
Indeed, already for intermediate doping, x = 1/9, 1/8 
or 1/5, spin-orbital order may change globally [12] ■ Here 
we investigate higher periodic doping x = 1/4, where 
half of the superexchange bonds are 3d-Ad hybrid ones 
and we show that they dominate and dictate the overall 
spin-orbital order. 

One finds that when at least one of the two impurity 
parameters is small, either Jj mp or ?7i mp (3), c doublon 
sites are doped and host-impurity coupling is AF. This 
AF coupling is amplified by charge excitations on 3d- 
Ad bonds when doublon orbitals at Ad ions are inactive 
it dominates in this range of the phase diagram and 
induces FM order along the Ad-Ad host bonds within 





Fig. 3 Ground state diagram obtained for periodic x = 1/4 
doping of d 4 spin-orbital host by d 3 ions (arrows). Dashed 
blue lines separate 2x2 unit cells in different phases; ques¬ 
tion marks in FSa2 indicate classically frustrated spins which 
order (small arrows) due to quantum fluctuations. 


the antiferrimagnetic (AFI) spin order for the entire 
ab plane, see Fig. 3. Thus only every second undoped 
vertical line is FM, as in the initial G-AF phase, while 
host spins are inverted on any other vertical line and 
the doublon orbitals flip from a to b. Indeed, this mod¬ 
ification of the orbital order stabilizes the FM inter¬ 
actions for c-a doublon pairs on the horizontal bonds. 
For FM order this orbital pattern is stabilized by dou¬ 
ble exchange, similar to a hole doped tig system [ 2 j. 
In contrast, if either Ji mp /Jhost or 7 ?i mp is bigger than 
~ 0.5, orbital dilution occurs on sublattice a. When the 
3d-Ad superexchange is AF, the orbital order along the 
undoped vertical FM lines is unchanged and one finds 
in phase AFa2 again the same spin order as in AFc, 
but now there are twice as many c doublons as a ones. 

As in the dilute limit (Fig. 2), impurity spins are 
frustrated also at x = 1/4 doping in phase FSa2 when 
?7i mp ~ ?7 ; c mp (Fig. 3) and the 3d-Ad superexchange nearly 
vanishes and changes sign. Then the host spins in be¬ 
tween experience almost entirely the AF Ad-Ad superex- 
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change and follow C-AF order. Quantum fluctuations 
stabilize opposite to them orientation of impurity spins, 
and this phase may be seen as a precursor of the FMa2 
phase which has again 75% FM bonds as the two AFI 
phases. At sufficiently large 77i mp FM order takes over 
in FMa2 phase (except for the range of c doublon dop¬ 
ing). The G -AO order is the same in this latter phase 
as in the undoped host, see Fig. 3. 


4 Discussion and summary 

Orbital dilution will play a role in several Mott insu¬ 
lators with spin-orbital order doped by ions without 
active orbitals. We have shown that the orbital order 
around impurities changes in general, so even in the 
dilute limit one may expect observable effects due to 
islands of reversed spins and doublon orbitals. 

The phase diagram of Fig. 3 confirms the general 
rule that only a sufficiently strong coupling Jrnip/^host 
leads to orbital dilution on a sublattice and to a rich 
competition between various types of spin order jT2j . 
We argue that the general trends found here are generic 
and the phase diagrams are only quantitatively modi¬ 
fied by quantum fluctuations, at least in systems with 
weak spin-orbit coupling. We have found that double 
exchange leads to local or global changes of spin-orbital 
order, similar to formation of orbital stripes [2] or or¬ 
bital polarons in doped manganites [4j. Such changes 
are expected to generate novel spin-orbital-charge mod¬ 
ulated patterns reported recently for f 2 g systems m- 

Summarizing, this study highlights the role played 
by orbital dilution due to d 3 impurities in cubic spin- 
orbital systems and opens a new avenue towards theo¬ 
retical understanding of Mn-doped layered ruthenates 
and related systems. We have shown that impurities 
change radically the spin-orbital order around them in 
the entire parameter range. As a general feature one 
finds frustrated impurity spins and their ability to po¬ 
larize the host orbitals around them. This property is 
remarkable and concerns both the dilute and high dop¬ 
ing regime, so one expects global changes of spin-orbital 
order in doped materials. It is challenging to investigate 
whether disordered impurities would generate similar 
changes of spin-orbital order as well. 

Acknowledgments W.B. acknowledges support by 
the European Unions Horizon 2020 research and inno¬ 
vation programme under the Marie Sklodowska-Curie 
grant agreement No. 655515. M.C. acknowledges fund¬ 
ing by the European Union FP7/2007-2013 programme, 
Grant Agreement No. 264098 - MAMA. We acknowl¬ 
edge support by Narodowe Centrum Nauki (NCN, Na¬ 
tional Science Center), Project 2012/04/A/ST3/00331. 


References 

1. Oles, A.M.: Fingerprints of spin-orbital entanglement in 
transition metal oxides. J. Phys.: Condens. Matter 24, 
313201 (2012) 

2. Wrobel, P., Oles, A.M.: Ferro-orbitally ordered stripes in 
systems with alternating orbital order. Phys. Rev. Lett. 
104, 206401(2010) 

3. Tokura, Y.: Critical features of colossal magnetoresistive 
manganites. Rep. Prog. Phys. 69, 797 (2006) 

4. Kilian, E., Khaliullin, G.: Orbital polarons in the metal- 
insulator transition of manganites. Phys. Rev. B 60, 
13458 (1999) 

5. Feiner, L.F., Oles, A.M.: Orbital liquid in ferromagnetic 
manganites: The orbital Hubbard model for e g electrons. 
Phys. Rev. B 71, 144422 (2005) 

6. Vernay, F., Penc, K., Fazekas, P., Mila, F.: Orbital de¬ 
generacy as a source of frustration in LiNi02- Phys. Rev. 
B 70, 014428 (2004) 

7. Reitsma, A., Feiner, L.F., Oles, A.M.: Orbital and spin 
physics in LiNKfU and NaNKU- New J. Phys. 7, 121 
(2005) 

8. Fujioka, J., Yasue, T., Miyasaka, S., Yamasaki, Y., 
Arima, T., Sagayama, H., Inami, T., Ishii, K., Tokura, 
Y.: Critical competition between two distinct orbital-spin 
ordered states in perovskite vanadates. Phys. Rev. B 82, 
144425 (2010) 

9. Khaliullin, G., Horsch, P., Oles, A.M.: Theory of optical 
spectral weights in Mott insulators with orbital degrees 
of freedom. Phys. Rev. B 70, 195103 (2004) 

10. Fujioka, J., Miyasaka, S., Tokura, Y.: Doping vari¬ 
ation of anisotropic charge and orbital dynamics in 
Yi-^Ca^VOa: Comparison with Lai-^Sr^VOs. Phys. 
Rev. B 77, 144402 (2008) 

11. Avella, A., Horsch, P., Oles, A.M.: Defect states and 
excitations in a Mott insulator with orbital degrees of 
freedom: Mott-Hubbard gap versus optical and transport 
gaps in doped systems. Phys. Rev. B 87, 045132 (2013) 

12. Brzezicki, W., Oles, A.M., Cuoco, M.: Spin-orbital order 
modified by orbital dilution in transition-metal oxides: 
From spin defects to frustrated spins polarizing host or¬ 
bitals. Phys. Rev. X 5, 011037 (2015) 

13. Qi, T.F., Korneta, O.B., Parkin, S., De Long, L.E., 
Schlottmann, P., Cao, G.: Negative volume ther¬ 
mal expansion via orbital and magnetic orders in 
Ca 2 Rui_ a ,Cr a; 04 (0 < x < 0.13). Phys. Rev. Lett. 105, 
177203 (2010) 

14. Mesa, D., Ye, F., Chi, S., Fernandez-Baca, J.A., Tian. 
W., Hu, B., Jin, R., Plummer, E.W., Zhang, J.: Single¬ 
bilayer E- type antiferromagnetism in Mn-substituted 
Sr3Ru 2 07: Neutron scattering study. Phys. Rev. B 85, 
180410(R) (2012). 

15. Cuoco, M., Forte, F., Noce, C.: Interplay of Coulomb 
interactions and c-axis octahedra distortions in single¬ 
layer ruthenates. Phys. Rev. B 74, 195124 (2006) 

16. Brzezicki, W., Dziarmaga, J., Oles, A.M.: Noncollinear 
magnetic order stabilized by entangled spin-orbital fluc¬ 
tuations. Phys. Rev. Lett. 109, 237201 (2012) 

17. Daghofer, M., Oles, A.M., von der Linden, W.: Orbital 
polarons versus itinerant e g electrons in doped mangan¬ 
ites. Phys. Rev. B 70, 184430 (2004) 

18. Brzezicki, W., Noce, C., Romano, A., Cuoco, M.: Zigzag 
and checkerboard magnetic patterns in orbitally direc¬ 
tional double-exchange systems. Phys. Rev. Lett. 114, 
247002 (2015) 



